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1 Introduction 



The interest in the existence and properties of nonhnear locahzed excitations 
(NLEs) in Hamiltonian lattices has considerably increased since the work 
of Takeno et al ([|ll,0 and references therein). The NLEs are classical so- 
lutions and represent localized vibrations of the constituents (particles) of 
the Hamiltonian lattice. NLEs are thus similar to breather solutions of the 
sine-Gordon equation (a partial differential equation) and sometimes called 
discrete breathers. Remarkably the lattice exhibits discrete translational in- 
variancy, i.e. no defects or disorder are present. Then the only reason for the 
possibility of the existence of NLEs is the nonlinearity of the system together 
with its discreteness. The nonlinearity is needed in order to be able to tune 
oscillatory frequencies by changing the energy or amplitude. The discrete- 
ness guarantees a finite band width of the frequency band of small-amplitude 
linearized extended excitations (phonons). Consequently one can avoide res- 
onances between the phonon band and the frequencies of NLE solutions. 

Recently the problem of finding time-periodic NLE solutions on Hamil- 
tonian lattices was reformulated into finding appropriate solutions of a set 
of coupled algebraic equations The variables of this set are the Fourier 
coefficients of the NLE solution. This approach opens several possibilities 
to obtain analytical and numerical results for NLEs. Especially spatial de- 
cay laws for every Fourier coefficient are derived, based on a linearization 
procedure 0. 

In the present work we will first formulate a new mapping in order to 
numerically solve the mentioned set of algebraic equations. We obtain high- 
precision numerical solutions for these equations. The level of accuracy al- 
lows us to check analytical predictions. Next we reconsider the spatial decay 
problem of the NLE solution and find, that under certain conditions the pre- 
viously applied linearization procedure does not hold without exceptions. We 
obtain nonlinear corrections and formulate conditions for their restricted ap- 
plicability. Finally we make use of the availability of high-precision numerical 
solutions and test our analytical predictions. 
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2 Formulation of the problem 

We consider a system of classical interacting particles arranged on a certain 
lattice. By that we mean that the positions of all particles in the groundstate 
(potential minimum) define a regular lattice structure. The potential energy 
of the system, expressed in the relative displacements of the particles from 
their groundstate positions exhibits a minimum if all displacements are zero. 
Since we are interested in the easiest demonstration of what we might call 
a general phenomenon, we restrict ourselfes to the case of a one- dimensional 
lattice with one degree of freedom per unit cell. Also we choose nearest 
neighbour interaction. We show later that these restrictions are at best of 
some technical use, i.e. more complicated systems will exhibit the same 
behaviour. The Hamilton function of our system can be then given in the 
form 



Here P and X denote canonically conjugated momentum and displacement 
of a particle and the label / counts the unit cells. The on-site potential V{z) 
and the interaction potential ^{z) can be choosen in a Taylor expansion 
representation: 




I 



(1) 




(2) 



(3) 



Finally the equations of motion are given by 

dH 



(4) 





(5) 




Fourier series representation for 



Xiity. 



+ 0O 



(6) 



fc=— oo 
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where uji = 27r/Ti. Let us separate the hnear terms in the displacements X 
on the right hand side of from the nonhnear ones: 



Xi = -V2X1 - 02(2X, - Xi_, - + Fr\X,) . (7) 

The nonhnear part F^*^"'^ of the force contains only higher than linear terms 
in the displacements: 



Fr=- E 

a=3,4. 



a-iy. ' {a-l)\ 

(8) 

Since all functions Xi{t) of our proposed solution (H) are periodic with the 
same period Ti, we can also expand F^*-"''' into a Fourier series: 

+00 

F/"')(t)= Fir'^e^'^-* . (9) 



The Fourier coefficients -F^"''' will be some functions of the coefficients Ak'i' 
from (depending on the given form of the potential functions (@),®). 
Inserting into (H) and using we finally obtain: 

^ujIAm = V2AM + M^^M - Ak,i-i - Ak^i+i) + Flf . (10) 

The infinite set of equations (|10|) (for all k) represents a coupled set of equa- 
tions for the Fourier coefficients A^i. 

This set of equations has to be analyzed in order to obtain nonlinear 
localized excitations, i.e. (cf. (^) solutions with 

/lfc,;^±oo^O . (11) 

Let us shortly review what is known about the solutions of (ITO) with condition 
(pHI). In it was shown, that the finding of a nonlinear localized excitation is 
equivalent to obtaining common points (homoclinic points) of two separatrix 
manifolds of a certain mapping. It was also shown there, that if a solution 
is can be found for a given system, then it is structurally stable with respect 
to perturbations of the Hamilton function and thus of generic type. Also 
in [H an existence proof of NLEs was given for a system with V{z) = 
and ^{z) = {l/2m)z'^^ (m = 2,3,...). Furthermore MacKay and Aubry 
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have shown, that in the hmit of weak interaction $(-2) and any anharmonic 
potential V{z) NLEs exist, remarkably for any lattice dimension, and also 
for larger interaction ranges and numbers of degrees of freedom per unit cell. 
Consequently the above considered simplest case is not different from more 
complicated (and therefore more realistic) models. A systematic study of 
one-dimensional 0,0 and two-dimensional lattices shows that there is 
absolutely no difference in the appearance and properties of NLEs. 



3 The method of finding solutions - designing 
a map 

Let us reformulate Eq. (^) into a map in order to obtain a numerical 
procedure for finding NLE solutions. For that we rewrite (^) in two ways: 



A 



kl 



1 



A 



kl 
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^k'l' 



V2 



{k 



2, ,2 



{i) 
kl- 



, Aii) > 
1 ^ ^k,l+li 



^{nl) 
kl 



{A 



k'V 



, (12) 
• (13) 



Here the notion of the upper labels {i + 1) and {i) can be ignored for a 
moment. For clarity we explicitely indicate in (|T2|) , ([T3|) the dependence of 



the coefficients -Ffc"'^ on the coefficients Aki. Both equations (|T2[),(|T3D are 



equivalent to each other with respect to solutions. If we now take into account 
the upper labels {i + 1) and (i), then these two equations in fact become 
two maps. Given the coefficients A^i at the z-th step, we can calculate the 
coefficients A^i at the {i + l)st step. If we find a fixed point of any of the two 
maps, then it will be also a fixed point of the second map. Each fixed point 
of these maps is a solution of the original equation ([l0|) . 

One fixed point can be immediately found - it is A^i = for all k,l. If we 
linearize the map around this fixed point and assume for a moment 02 = 
we obtain the eigenvalues 



A 



(12) 
kl 



V2 



A 



k^Ujf 
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V2 A® 



(14) 
(15) 
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Here the upper labels of the eigenvalues indicate the equation number of the 
map they belong to. If now 02 is smoothly increased from zero, then the 
eigenvalues (^l|),(|I3p will also smoothly change 0. Consequently we might 
expect some control over the values of the eigenvalues in the case of nonzero 
interaction. As for the eigenvectors, any orthonormal set of basis vectors 
would be an eigenvector base for the linearized maps in the case 02 = 0. For 
clarity we can always choose the eigenvector set which appears in the limit 
02 0. Then the eigenvectors are also smoothly changed by increasing 02 
from zero. 

As it follows from both eigenvalues for a given pair of {k,l) are 



positive and inverse to each other. Consequently one of the eigenvalues will 
be larger than one and the second eigenvalue will be smaller than one. In 
order to find a NLE solution we design a new map out of the two maps 
(|T2p,(|T3p. This new map can be composed out of the two maps (P!^),(P^) 
in an appropriate way depending on the given problem one wants to solve. 
Let us explain the strategy of this designing. Suppose we want to find a 
NLE solution with a given frequency cui, which is centered on a lattice site 
/ = 0. Then we choose out of the two maps (|l3)?(IIlD the one which yields an 
eigenvalue larger than 1 for the pair {k = ±1,1 = 0) according to ([T^) , ([T3|) . 
For all other pairs {k, I) we choose the map which yields an eigenvalue smaller 
than 1 (note that this choice might be still different for different k). By that 
we achieve the following. If we choose as an initial condition for our designed 
map a set of the A^i with Ak=±i,i=o = S, Akj^±i,i^o = 0, then we can expect 
that starting iterating our map we essentially get growth in the direction of 
the unstable eigenmode at {k = ±1,1 = 0), whereas all other eingenmodes 
are not excited (eigenvalues lower than one). Because we have designed a 
map in which the translational invariancy with respect to I is broken at Z = 0, 
we also know that the corresponding eigenvector is of local character for small 
02. A way to estimate whether this statement is correct is to calculate the 
ratio of the off-diagonal terms in the linearized map (which are roughly given 
by 02/^2) over the difference between the diagonal terms: 

^ _ ^ ^ ^ ^ (16) 

f 2 iOi f 2 f 2 — t^l 

This ratio has to be small enough in order to be sure that the above given 
arguments hold. Once our iteration yields local growth for the Fourier co- 
efficients, we can hope that the nonlinear terms of the map, which become 
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important for large enough A^i, will lead the iteration to the fixed point 
which corresponds to an NLE solution. There is absolutely no guarantee for 
that to happen. We have no knowledge about the properties of this fixed 
point, even if we would know for sure that the NLE solution exists (which 
means in turn that the fixed point exists). All the designing is focussing on 
is to get the right initial growth (namely local growth) in our iteration. 

If we would search for a NLE which is centered between two lattice sites 
I = 0, 1, all we would have to change is to add an unstable map for the 
pair k = ±1,1 = 1, and depending on the symmetry of the solution (in- 
phase or out-of-phase motion) to choose equal or opposite signs for the initial 
perturbations in Ak=±i,i=o and Ak=±i,i=i. 

4 Solutions 

We implemented two numerical realizations for the designed map. The main 
difficulty is to account for the nonlinear terms -F^"''' (^). In this paper we 
report on results for systems with V{z) and ^{z) being finite order polyno- 
mials in z. Then one way of calculating the nonlinear terms is to evaluate all 
possible combinations of products of the Fourier coefficients Aki. We will call 
this method polynomial approach. As an example we give the contribution 
of the first term in the sum on the right-hand side of equation (|]): 

^kl = X! / _ XI Akj^iAk2l...Ai:^^-^iSkXki+k2+...+kc-i) ■ 

o=3,4,... v /■ ki,k2,...,ka—i=—oo 

(17) 

Here Sa,b stand for the Kronecker symbol with integer a, b. 

The polynomial approach can become very unefficient if the number and 
order of nonlinear terms in ([T7|) become large (in fact a = 5 can already 
pose serious computing problems). If the potential functions are not given 
by finite order polynomials then the polynomial approach breaks down. In 
this case it is more useful to numerically integrate the function Fi^^\t): 

Flf = 1 r\ F/"')(t)e-^^-Mt . (18) 

il J-T/2 

We will call this approach Fourier integral approach. Here Ti = 27r/ci;i is 
the period of the NLE solution. If we know A^i for all {k,l), then we can 
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immediately calculate ( ]I8| ) using (||). In order to test the accuracy of our 
calculations we can use both approaches and thus compare the two results, 
which have to be equivalent. 

Another important restriction on the numerical approach is the finiteness 
of the /c-sums. We have to introduce by hands a cutoff in /c-space: A^i = 
for \k\ > kjnax- This cutoff will be justified only if the numerically evaluated 
Fourier coefficients A^i will drop to small enough values at the cutoff k^ax- 
Here small enough will be defined through the accuracy of the numerical 
calculations. We will stop the iteration if the sum 

El4^-4r'^l<io-^° ■ (19) 

k,l 

In this paper we will study three models. Model I is characterized by 
V{z) = 1/2^2 _ 1/4^4 g^j^^ ^(^^) ^ l/2Cz\ Model II has the potential 
functions V{z) = - + 1/Az\ ^{z) = l/2Cz^. Model III is defined 
through V{z) = l/2z^ + l/Az^, ^{z) = l/2Cz^ . Model I is characterized 
by a local minimum for z = and a global instability for \z\ > 1. Here 
we will restrict ourselfes to small enough amplitudes such that we avoid the 
instability region. Model II is the well known <l>^-lattice model (after a change 
of variables z' = z — 1 one arrives at the symmetric form of the double-well 
onsite potential V{z') = l/4{z''^ — 1)^). This model has in fact two identical 
groundstates z' = ±1 or z = 0, 2. Without loss of generality we have choosen 
2;' = — 1 as the groundstate around which we expand the potential functions. 
Model HI differs from Model I by the sign in front of the quartic term in 
V{z). Consequently there is no instability in this model. In all models the 
parameter C = 02 regulates the bandwidth of phonons obtained after a 
linearization of the equations of motion around the groundstate: 

V2 < ojI <V2 + AC . (20) 

Here Uq is the frequency of a phonon wave with wave number q. For our 
models we have: model I - f2 = 1; model II - f2 = 2; model HI - f2 = 1. In 
our calculations we will use C = 0.1. 

In all three model cases the phonon band has a nonzero lower band edge 
f2. In order to solve the map for the Fourier coefficients we have to know (or 
be able to predict) whether the NLE frequency can be below or above the 
phonon band (or may be below and above). There exists a rather straightfor- 
ward method to do so. First we have to know the symmetry of the solutions 
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with respect to their position on the lattice. In most cases (even supported by 
rigorous proofs for some cases) we can expect two solutions: a NLE centered 
on a lattice site (meaning that the center of energy of the energy distribution 
of the solution has a maximum on a lattice site) which we call XC solution; 
and a NLE centered between two lattice sites which we call XN solution, 
(this notation was first used in |jlO[)- Now all one has to do is to consider two 
possibilities: that the NLE frequency might be below or above the phonon 
band. As it was shown in P], this option defines whether the motion of 
particles in the solution is of in- or out-of-phase character. In 
explained how to construct a certain effective potential out of these input 
informations. The motion of a particle in this effective potential can be then 
analyzed with respect to the energy dependence of its oscillation frequency. 
This energy dependence is qualitatively (and often even quantitatively) the 
same as the energy dependence of the NLE. If e.g. the initial assumption 
was 'frequency below the phonon band', but the effective potential yields 
only frequencies above the phonon band, then the initial assumption was 
wrong (and vice versa). Let us give the form of the effective potentials for 
the 4 different cases: 

V,ff{z) = V{z) + ^z) + ^-z) , XC , uoi<u, ,(21) 

V,ff{z) = 2V{z) + ^{z) + $(-z) , XN , iJi<iUg , (22) 

Veffiz) = V{z) + ^z) + <l>(-^) , XC , , cui > cu, , (23) 

Veff{z) = V{z) + V{-z) + ^2z) + ^z) + ^-z) , XN , ooi > .(24) 

As it follows from the analysis, model I allows only for cases (pll),(|22|) {ui < 
ujq). For model III we find the reverse case uji > Uq. For model II and small 
energies we find uJi < Uq, whereas for large energies also ui > ujq can be 
realized. 

Having this information, we can proceed with the numerical treatement of 
the map for the Fourier coefficients. In Fig.l we show an XC NLE solution as 
obtained for model I and the frequency uji = 0.8. We used periodic boundary 
conditions (PBC), 100 lattice sites and kmax = 30. Because of the symmetry 
of the potential functions all Fourier coefficients with even k vanish. In our 
calculations (both with the poljTiomial and Fourier integral approach) we 
however did not use this information for the calculation. Both methods yield 
exactly the same solution, in particular the even Fourier coefficients indeed 
vanish. In Fig.l the dependence of \Aki\ on the lattice site / is shown in a 
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semi-logarithmic plot. The actiial data are given as open squares. Data for 
equal values of k are connected with solid lines. The value of k increases from 
top to bottom: A; = 1, 3, 5, .... As is clearly seen in Fig.l, Fourier coefficients 
for each k show up with a logarithmic decay in /. What is not observable in 
the plot of Fig.l is that the Fourier cocfficcnts for a given k alternate their 
signs as one increases or decreases / by 1. Only for k = 1 the sign of all 
All is the same, independent on /. The Fourier coefficients for negative k 
are identical with the corresponding coefficents for positive k: A^i — A_k,i- 
Let us also mention, that we have checked any finite size effects. We have 
computed the same solution for 30 lattice sites, and found practically no 
difference in the values of the Fourier coefficients. This appears to be logical, 
since the solution we find is very localized, so that it does not matter, how 
large the system is. Consequently we can be sure that we find a solution 
which is also a solution of the infinite problem. 

In Fig. 2 we show an XN NLE solution for the same model I, where all 
other parameters are as in case of the solution in Fig.l. All the comments 
from Fig.l apply to Fig. 2. 

Let us consider model IL We choose a frequency a;i = 1.3 which is below 
the phonon band. Because of the assymetry of the corresponding onsite 
potential V{z) Fourier components with even k will be nonzero. The result 
for an XC NLE is shown in Fig.3. Data belonging to one /c-value are again 
connected with solid lines. The value for k from top to bottom is k = 
1, 0, 2, 3, 4, .... The signs of Aki are positive for A; = 0, 1, 3, negative for A; = 2 
and alternating for all other /c-values. We also created a similar XN NLE 
solution with same frequency for model II which behaves similar to the XC 
solution in Fig.3. 

In Fig. 4 we show the result for an XC NLE solution for model III. Here 
the frequency cui — 1.3 is choosen to be above the phonon band. As in the 
case of model I the symmetry of the potential functions causes all even k 
Fourier coefficients to vanish. The fc-values for the nonzero odd k Fourier 
coefficients in Fig. 4 increase from top to bottom : k = 1,3,5,.... For all 
odd A;-values the signs of the Fourier coefficients alternate as we change the 
lattice site Z by 1. Again we also created a similar XN NLE solution. 
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5 The spatial decay revisited 



In [0] the spatial decay of the Fourier components Aki was studied. The main 
idea was to start with (p!oD and rewrite it in the form 

Ak,i+i = ^Aki + 2AM-Ak.i-i-^Aki + ^Fl:f . (25) 

02 02 02 

Equation (|25|) can be interpreted as a map: knowing all Fourier coefficients at 
lattice sites I and (^ — 1) one can calculate the coefficients at + Naturally 
one can iterate also to negative values of / (i.e. one just exchanges (/ + 1) with 
(/ — 1) in (^)). If we iterate into the tails of a NLE solution, then because 
of ( PTD it was suggested in that the nonlinear contributions -F^'^''' which 
are at least of second order in Ak'u can be neglected. Then equation (|^) is 
transformed into an infinite set of two-dimensional maps for each value of k. 
As a result exponential decay laws are obtained 0: 

A,, ~ (sgn(A,))' e^^^l' , A, = 1 + ^ ± ^(1 + ^)^-1 (26) 

with Kk{uJi) = (f2 — kluf) / {(f)2) ■ The sign in (|26|) hase to be choosen in order 
to obey the condition |Afe| < 1. One of the interesting properties of (p^) 
is that all one has to know about the problem is the frequency of the NLE 
solution and the positions of the phonon band which are given by (^2,02)- 
Consequently several predictions were drawn in on this basis which were 
successfully checked. 

However in order to justify the dropping of the nonlinear contributions 
in (^) one has to insert the predicted decay laws as obtained from the 
linearization into the full equations and show that the nonlinear terms are 
indeed negligible. Let us consider the contributions of the nonlinear terms as 
given in (0). First we note that the function d{k) = ln(|Afc|) is by definition 
always negative. d{k) is defined on a discrete lattice (integer k) and d{k) = 
d{—k). Let us consider k > 0. Then d{k) has a single maximum either for 
one value of k = km or for two subsequent values k = km and k = {km + 1). 
For k > km the function d{k) is monotonically decreasing with increasing 
k. Then the linearization in (|25|) is violated for a given ko if there exists a 
sequence of ki, k2, ka-i such that 

d{ko) < {d{ki) + d{k2) + ... + d{ka-i)) (27) 
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holds together with the condition (see ([T7|)) 

ko = ±ki±k2±...±ka-i . (28) 

The reason is that in this case a given nonhnear combination is decaying in 
space weaker than the corresponding hnear term. 

Since equations (p7|) ,(^8D can never be satisfied for ko = km it follows 
that the linearization is always justified for km- Let us consider k > km- If 
the sequence of {ki, ^2, ka-i) yields violation of the linearization, then all 
fc-numbers from this sequence have to be of lower value than k. Increasing k 
step by step starting from km we can then account for all possible nonlinear 
corrections of the result of linearization. Since d{k) ~ —2ln{k) for large 
values of fc, it follows from a numerical analysis that for each value of a there 
will be a fc+ such that no nonlinear corrections appear for /c > Also for 
a given function g{k) there will be a a+ such that no nonlinear corrections 
appear for any value of k. 

Let us consider the numerical results for NLEs from the previous section 
in the light of the above given considerations. The case shown in Fig.l 
was analyzed in 0. Fig. 5 in P shows very good agreement between the 
fc-dependent exponents of the decay as found from the numerical solution 
(Fig.l) and the result of the linearization. A careful check of the nonlinear 
corrections indeed show that the linearization result holds without exception. 
The same result is valid for the XN NLE in Fig.2. 

The XC NLE solution shown in Fig. 3 can not be described in its decay 
properties by the linearization result only. In Fig. 5 we show the exponents 
as calculated from Fig. 3 (squares) and the predicted exponents from the 
linearization procedure (filled circles). Clearly for k = 2 (and also less visible 
for k = 3) the numbers differ. The exponents for the solution are 



k num. result linearization 






-1.3202 


-1.3415 


1 


-0.6904 


-0.6898 


2 


-1.3796 


-1.6588 


3 


-2.0748 


-2.1143 


4 


-2.3957 


-2.3951 


5 


-2.6018 


-2.6026 


6 


-2.7663 


-2.7682 
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Nonlinear corrections can appear from terms with a = 3,4. For a = 3 we 
can find the following corrections. For k = 2 the nonlinear contribution 
^1/^1/ yields an exponent (—0.6904 — 0.6904) = —1.3796, which is exactly 
what is found from the numerical result and is larger than the value of the 
linearization result —1.6588. The sign of A21 should alternate according to 
the linearization result. If the nonlinear terms in (pSj) dominate, the sign of 



A21 ~ —All should be negative - exactly as found from the numerical result 
(cf. previous section). For k = 3 the nonlinear contribution A11A21 (with 
the corrected exponent for A21) yields an exponent (—0.6904 — 1.3796) = 
—2.0700 which is again very close to the numerical result and larger than 
the linearization result —2.1143. The sign of A31 should alternate according 
to the linearization result. If the nonlinear terms in (ESI) dominate the sign 
of ~ —A11A21 should be positive {An > is found from the numerical 
solution) - exactly as found from the numerical result (cf. previous section). 
The only nonlinear contribution from the nonlinear term with a = 4 appears 
for = 3 (^f;) and is equivalent to the a = 3 contribution {A11A21 = AuAfi) 
with respect to the value of the exponent as well as with respect to the sign of 
the correction. For > 3 no nonlinear corrections are found - consequently 
the linearization result coincides with the numerical findings both in the value 
of the exponent and in the sign of the Fourier components (alternating). 

The XC NLE solution shown in Fig. 4 has also nonlinear corrections in 
the spatial decay. In Fig. 6 we show the exponents as calculated from Fig. 4 
(squares) and the predicted exponents from the linearization procedure (filled 
circles). For k = 3 the results are different. The exponents for the solution 
are 

k num.result linearization 
1 -0.6722 -0.6709 
3 -1.9910 -2.1464 
5 -2.6103 -2.6133 
7 -2.9114 -2.9117 
9 -3.1324 -3.1325 



Only nonlinear terms with a = 4 are present. For k = 3 the nonlinear 
contribution AuAuAu yields an exponent (-0.6722 - 0.6722 - 0.6722) = 
—2.0166 - very close to the numerical result and larger than the value of the 
linearization result 2.1464. Since the sign of An alternates, so does the sign 
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of A^i ~ All. For all other k ^ 3 the nonlinear terms give no corrections to 
the linearization result. 

Nonlinear corrections can be expected if one of the multiples of the fre- 
quency kuji comes close to the phonon band - then the corresponding spatial 
decay exponent comes close to zero and can easily contribute to the decay 
properties of Fourier coefficients with e.g. k = {a — l)k. 



6 The decay in /c-space 

So far we have considered the spatial decay of NLE solutions. However as 
it was already explained, in order to obtain numerical solutions we have to 
introduce a cutoff in fc-space. This cutoff will be only justified if the Fourier 
components decay in the fc-space. From the examples in Figs. 1-4 this decay 
is clearly visible. In this section we will obtain some analytical results for 
the decay in /c-space for a solution where no nonlinear corrections to the 
linearization result (of the spatial decay, see previous section) appear. 

As it appears from (^) the fc-dependence of the Fourier coefficients is 
given by 

Am ~ s(fc)e-i'^(l^*l)l'l (29) 

where the function s{k) describes the fc-dependence of the Fourier coefficients 
at lattice site / = (for convenience we assume that the NLE is created with 
center at / = 0; in case of reference to numerical results a shift in the lattice 
site numbers is implied). If we consider the XC NLE solution in Fig.l, then 
we observe that ln{s{k)) ~ —k. For large values of k 



X 

and consequently 



\Aki\ ~ k~^^'^s{k) . (30) 

In case of Fig.l this implies that the decay of the Fourier components in 
/c-space is given by a product of an exponential decay and a power law decay, 
where the power law decay is lattice site dependent - the farer from the 
center of the NLE, the larger the exponent of the power law decay. In fact 



according to (|30D the exponent of the power law decay should increase by 2 



as one switches from a lattice site to its nearest neighbour. 
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In Fig. 7 we plot the normalized Fourier components \Aki\/ s{k) from Fig.l 
as a function of in a log-log plot. Data with same values of / are connected 
with lines. Clearly power laws are observed. In Fig. 8 the /-dependence of the 
power law exponent is shown by plotting the differences between power law 
exponents on neighbouring lattice sites - as predicted this difference appears 
to be approximately —2 (deviations for larger values of I appear because of 
the smallness of the unrenormalized Fourier coefficients). 

From the results in this section we can follow that the decay of the Fourier 
components in fc-space is as faster as the distance from the NLE center in- 
creases. 

7 A variational approach 

So far we have implemented a numerical method of obtaining NLE solutions 
and have used the numerical results in order to verify and improve several 
analytical findings. For completeness we show in this section that NLE so- 
lutions can be obtained out of a variational problem, i.e. we can define a 
function of the Fourier coefficients such that a NLE solution will correspond 
to an extremum of this function. This result can be useful both in consid- 
ering an independent method of finding numerical NLE solutions as well as 
in performing existence proofs for NLEs in lattices with dimensions higher 
than one and in considering the analogous properties of quantum lattices. 
Let us consider the action 5* of a given classical lattice: 



and the Hamilton function H is given in (|I]). We search for a periodic solu- 
tion (closed trajectory) with period t = iT^juox. Inserting (^ into ( PT| ) and 
integrating over the period we find for the action of a periodic solution 




(31) 



where the Lagrange function L is given by 



(32) 



I k 



oo 



1 



EE 




(33) 



I a=2 



kl,k2,...,ka 
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Here the ... in (^) implies the evaluation of the interaction terms which 
look similarly to the term in the second row of Eq. (^) except that the 
lattice site labels in one product can come from two neighbouring lattice 
sites. 

Let us consider the derivative of Sp with respect to a given A^i'. 

TT^^-T" = ^2t^M-fc,/-X! """T TYT X! ^kiiAk2i---Ak^_^iSo,{k+ki+k2+...+k^ 

^^^Am - V- k^,k2,-,ko,-i 

(34) 

Now we observe that if a periodic solution on the lattice is found, i.e. if 
equation (|^) holds, then for these solutions the right hand side of (|^) 
vanishes (if one replaces by — in {^)). Consequently the function Sp 
from (|33| ) has to have an extremum with respect to variations of the Fourier 
coefficients A^i for a periodic solution. Also the inverse is true: if Sp has an 
extremum with respect to variations of the Fourier coefficients 

- (35) 



dA 



kl 



for all (kl) then the set of Fourier coefficients which corresponds to this 
extremum is a periodic solution of the equations of motion. 



8 Discussion 

We have shown the construction of a map for the Fourier coefficients of a peri- 
odic solution on a lattice. This map is specifically designed in order to obtain 
spatially localized time-periodic solutions. The usefulness of this method is 
first that very little (in fact nothing specific) has to be known about the so- 
lution one searches for, i.e. the initial condition (starting distribution of the 
Fourier coefficients) has to reflect only the localized character of the solution 
one looks for. This circumstance makes our approach much more easy than 
e.g. the approach used in |T^. There also a map was designed in order to 



obtain localized time-periodic solutions. However the initial condition had 
to be choosen to be very close to the flnal result. In other words one had to 
obtain information about the flnal solution before even starting the map. 

Secondly with our method we derive flnal solutions with high numer- 
ical precision. This precision allows for the numerical analysis of several 
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problems, e.g. the dynamical stability of the derived solutions, the study of 
scattering events of plane waves by the localized solutions and others. Let us 
explain this using the scattering problem. In |TT| the scattering (reflection) of 
phonon waves by a single localized solution was studied for a one-dimensional 
lattice. The maximum amplitude for the localized solution was of the order 
of 1, and the errors of the amplitudes for the same solution were of the order 
of 10~^ leading to uncertainties in the energy distribution of order 10~^. The 
infalling phonon wave had to have an energy density of the order of 10~^ 
in order to ensure nearly linear dispersion. The transmittion of phonons 
was increasingly suppressed with increasing wave number. Since the squared 
transmittion coefficient was observable only for values down to 10~^ — 10~^ 
the maximum wave number under study was around OAqzB (here qzB means 
the zone boundary wave number). With the accuracy of the final solutions 
as derived in this work (10~^^ in amplitudes, 10"^'' in the energy density) one 
should be able to trace the squared transmittion coefficient down to values 
of 10~^° — 10~^^ which implies that nearly the whole Brillouin zone can be 
studied. 

Of course the presented method can be used in lattices with dimension 
larger than one too. The analysis of the solutions however will be much 
more tedious, thus we have demonstrated the method using one- dimensional 
systems. The formulation of the variational approach in the previous chapter 
can be also useful in lattice dimensions higher than one. 

We have used the results obtained with the new method in order to test 
predictions of analytical considerations of the spatial decay of a localized 
vibration. As a result we found that nonlinear corrections of a lineariza- 
tion method used to account for the exponential decay appear if a multiple 
of the fundamental frequency of the localized solution comes close to the 
phonon band. We have further shown how these nonlinear corrections can 
be obtained analytically and found very good agreement with the numerical 
results. This part of the work was also carried out for one-dimensional sys- 
tems. It is possible to analyze higher dimensional lattices with this method 
too. For that one has to consider lattice Green's functions (see e.g. [p!^,||TB|). 
These functions will appear if one deals with the equations for the Fourier 
components of a localized vibration on a given lattice (the equations are es- 
sentially as in (p!0|)). Far away from the center of the NLE one can again 
linearize these equations. The lattice Green's functions can be used then in 
order to obtain analytical solutions. The fc-dependent spatial decay will be 
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obtained, and nonlinear corrections can be considered in full equivalence to 
the cases analyzed in the present work. 

Let us finally stress that the method and results obtained in this work 
can be used in order to perform a dynamical stability analysis of the NLE 
solutions as indicated in The corresponding eigenvalue problem can be 
formulated in the space of Fourier coefficients. As an input of this stability 
analysis one needs precise data on the Fourier coefficients of the periodic 
NLE solution, which are provided by the presented method. 
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FIGURE CAPTIONS 



Fig.l 

Numerical result for an XC NLE solution for model I (see text). 
Open squares - numerical data. 

Lines connect Fourier coefficients with same value of k. 



Fig.2 

Same as in Fig.l but for an XN NLE solution. 
Fig.3 

Same as in Figs. 1,2 but for an XC NLE solution of model II. 
Fig.4 

Same as in Figs. 1,2, 3 but for an XC NLE solution of model III. 
Fig.5 

Exponents of the spatial decay ln|A| as a function of k for the NLE solu- 
tion shown in Fig.3. 
Open squares - numerical result. 
Filled circles - linearization result. 

Fig.6 

Same as in Fig.5 but for the NLE solution shown in Fig.4. 
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Fig.7 



fc-dcpendence of the Fourier coefficients of the NLE solution shown in Fig.l. 
Here the Fourier coefficients are renormahzed using the data for the central 
particle. The result is plotted in a log-log plot. 
Open squares - numerical data. 

Lines - best eye fit of a power law (straight line in the plot) for every lattice 
site I. 

Fig.8 

The difference of the power law exponents for the /^-dependence of the Fourier 
coefficients for neighbouring lattice sites as obtained from Fig.7. 
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